Fractal fluctuations in quantum integrable scattering 
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We theoretically and numerically demonstrate that completely integrable scattering processes 
may exhibit fractal transmission fluctuations, due to typical spectral properties of integrable sys- 
tems. Similar properties also occur with scattering processes in the presence of strong dynamical 
localization, thus explaining recent numerical observations of fractality in the latter class of systems. 
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Scattering processes dominated by a statistically large 
number of metastable states display reaction rates which 
depend on energy or other control parameters in com- 
plicated ways. Reaction curves (i.e., curves obtained 
by plotting cross sections versus the relevant parame- 
ters) demand statistical description. This approach, orig- 
inated in nuclear physics, has given rise to the theory of 
quantum chaotic scattering 0] . Classical chaotic dynam- 
ics plays an important role therein [0 , with applications 
in many fields, including the study of mesoscopic conduc- 
tance fluctuations |l . 

In the presence of complicated reaction curves, a natu- 
ral question is about their fractality. The latter is not 
meant in the strict mathematical sense, but rather as 
a property to be possibly observed on a wide range of 
resolution scales. Reaction curves probe the location of 
resonance poles in the complex energy plane, so their 
fluctuation properties - including fractality - are encoded 
in the distribution of such poles. In particular, the lack 
of smoothness required by fractality can only be pro- 
duced by resonance poles clustering in the vicinity of the 
real energy axis. This excludes fractal reaction curves 
in quasi-classical cases when the underlying classical dy- 
namics is completely chaotic (i.e, uniformly hyperbolic), 
because the quantum resonances are then concentrated 
away from the real axis. This reflects the exponential 
decay in time exhibited by such systems over long time 
scales. A dynamical signature of poles clustering near the 
real axis is instead slow algebraic decay of the survival 
probability inside the interaction region. On the classi- 
cal level, such slow decay is in particular exhibited by 
systems with a mixed phase space, endowed with a hier- 
archical structure of stable islands 0. On such grounds, 
fractal fluctuations were predicted by Ketzmerick for the 
quasi-classical transmission of electrons through meso- 
scopic cavities j5|. This prediction has received numer- 
ical I! and experimental |7| support. However, fractal 
fluctuations have been numerically observed also in 2d 
tight-binding models of quantum dots, where the relation 
to classical dynamics is unclear H, and even in models 



where this relation is irrelevant ||, due to strong quan- 
tum localization on one hand, and to absence of signifi- 
cant classical critical structures on the other. It therefore 
appears that fractality of reaction curves is not strictly 
associated with critical structures in the classical phase 
space. 

The issue of fractality in the remaining major class of 
dynamical systems, namely the integrable ones, has not 
yet been investigated |[(|. This is the main purpose of 
this Letter. We first identify special statistical properties 
of resonance poles, which afford fractal reaction curves. 
This we do in absolute generality, without assuming inte- 
grability or other special properties of the quantum sys- 
tem. We then argue that such properties are optimally 
exhibited by suitable, completely integrable processes, 
thanks to a well known generic property of integrable 
systems, that their energy levels share some of the prop- 
erties of a random sequence [pd| . We use a textbook 
scattering problem to confirm our theory with numerical 
data. Finally we note that the same properties are also 
typical of fully chaotic, yet strongly localized systems, 
and thus we explain the results of ref. |J. 

We consider a weakly open quantum system, with scat- 
tering resonances at Ej — iTj/2 and the real energies Ej 
arranged in increasing order. The energy dependence of 
a typical cross section T(E) consists of a smooth back- 
ground plus a resonant part, which we write in the form: 
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with Cj slowly varying with j. We restrict within an en- 
ergy interval (E — W, E + W), and we assume e -C W <C 
Eq, where e is the mean level spacing. To perform fractal 
analysis of the graph of T r (E) vs E we divide the interval 
(E -W,E + W) into subintervals A fc) (k = 1, 2, M) 
of equal size 8 cx 1/M. Upon each subinterval we pile up 
squares of side 6, and denote by N(S) the total number 
of squares met by the graph of T r (E) . Algebraic scaling 
N{8) cx with / > 1 between scales 5 m i n <C S max 
signals that in between such scales the graph exhibits a 



1 



fractal (box-counting) dimension /. In order to deter- 
mine / we compute 

\og(N(S)) _ \ogjM- 1 J2 k a k (T r j) _ \og{a(T r )) 
log(^i) ~ + log(M) ~ log(<5) 

where <Tk{T r ) > is the maximal excursion of T r {E) 
within the k— th interval. We next assume the follow- 
ing. First, in the given energy interval, the frequency 
P(T) of widths less than T scales like P(T) ~ aT 1 ^ 01 
at e < r < r, with T the mean resonance width, and 
< a < 1. Second, both Tj and Ej form uncorre- 
lated sequences. Finally, resonances are strongly over- 
lapped, in the sense that e -C T <C W. Then, at 
e <C 5 <C r, the fluctuation of T r in an interval A& is 
mostly due to many tiny resonant peaks which are cen- 
tered inside the interval and are narrower than S. Every 
such peak contributes a Lorentzian term in ([!]), and the 
mean square oscillation of this term as E ranges in is 
~ S^PiSy 1 f* dP(T)T ~ 1. There are n{6) ~ e^SPiS) 
such peaks; as they contribute uncorrelated oscillations, 
we estimate (a(T r )) ~ y/n(S) ~ e^S 1 -"/ 2 . There- 
fore, / = 1 + a/2 in a range S min < S < 5 max , with 
S m ax ~ r, and Smin roughly estimated by n(S min ) ~ 1, 
that is, S mm ~ (e/a) 1 /' 2 ""). 

Fluctuating cross sections may also be generated at 
fixed energy, by varying other parameters, as in the case 
of magnetoresistance fluctuations in mesoscopic physics. 
The above analysis carries over to such fluctuations on re- 
placing energy by the relevant parameter, provided that 
the above assumptions remain satisfied. 

The above described conditions are met in some phys- 
ically relevant situations. Uncorrelated sequences of en- 
ergy levels are a distinctive feature of generic integrable 
systems Jll]] . We then surmise fractal reaction curves for 
completely integrable scattering processes, provided they 
display a slow decay, leading to an inverse power law dis- 
tribution of T's. We shall presently describe an explicit 
example of such a process. 

We consider the quantum dynamics of a particle of unit 
mass moving inside the infinite strip < x < L x in the 
{x,y) plane, with tunneling barriers at y — ±L y /2. The 
Hamiltonian is 

h 2 

H = -yA + had(y - L y /2) + haS{y + L y /2) 

We use periodic boundary conditions at x = 0, x = L x . 
The Dirac delta functions enforce additional boundary 
conditions: d y ip(x,y+) — d y tp(x,y—) = 2ah~ 1 'ip(x,y) at 
y = ±L y /2 . The physical model is a rectangular bil- 
liard, whence the particle can escape into semi-infinite 
leads, by tunneling through the two horizontal sides. For 
\y\ > L y /2 eigenfunctions are u% m {x,y) = m (x)Of(y), 
<j> m (x) = exp(2nimx/L x ), 0£(y) = A ± (k) exp(ik\y\) + 
B ± (k)exp(-ik\y\) and 2E/h 2 = k 2 + Ait 2 m 2 / L 2 X , with ± 
denoting the upper (y > L y ) and the lower {y < —L y ) 



lead respectively. For given E > there are a finite 
number of open scattering channels labelled by the in- 
teger m, \m\ < lnt{L x {-nTi)~ l \J E /2) , and by the lead 
label ±. The coefficients A il {k) are related to B ± (k) by 
the scattering matrix. Scattering resonances are located 
at complex values of energy z„ jm = E niin — iT n /2 = 
2i: 2 m 2 h 2 /L 2 + k 2 h 2 /2, (n = 1,2,...), where k n arc the 
complex roots of the equations: 

e ikLy ± J -p ikh( j-l = (2) 

The following asymptotic formulae can be computed by 
using Lagrange's theorem on the inversion of analytic 
functions: 

e~.„ = <™ - § { t arcta °<^ » + (A) } m 

where v n = nirh/L y is the velocity in the n— th vertical 
mode of the closed rectangle, and En) m = 2ti 2 -K 2 m 2 / L 2 + 
v 2 /2 are the eigenvalues of the closed (<r = oo) billiard. 
The 1st term on the rhs of the 2nd eq.(^) is the decay 
rate of a classical billiard ball inside the closed rectan- 
gle, with velocity v n in the y— direction, and absorption 
probability at y = ±L y /2 equal to the transmission coef- 
ficient for a plane wave through on a 6 barrier, given by 
v n( a2 ~ >rV n)~ 1 - We assume h <C oL y and thereby neglect 
0(fi/aL y ) corrections in (^). The statistics of real parts 
of resonances is then like the energy level statistics of an 
integrable system, with a mean level spacing only dif- 
ferent by corrections of order (fr/<jL y ) 3 from that of the 
closed billiard: e = 2irft 2 {L x L y )~ 1 . The correction to the 
closed billiard levels in (||) lifts possible degeneracies due 
to commensurate geometry. 

In an energy interval (E - W, E + W), e < W < E , 
the smooth (Thomas-Fermi) part of the integrated P(T) 
distribution is computed from a microcanonical distribu- 
tion of classical billiard trajectories at energy E , each 
with a decay rate T (||). This quasi- classical distribu- 
tion has a mean r^; ~ ha~ 2 L y 1 E^ 2 , and behaves like 
aT 1 / 3 at small T, with a = (vr 2 E /2)- 1 / 2 (L y a 2 /K) 1 / 3 . 
The same behaviour can be assumed for the quantal dis- 
tribution, provided that Eq < a 2 , and that e <C r^o- 
The former condition ensures T 1 / 3 behaviour of the quasi- 
classical distribution in a range of T— values comparable 
to Te. Together with the latter condition it ensures that 
the quasi-classical behaviour is observed over a statisti- 
cally significant number of resonances. Reflecting generic 
properties of energy spectra of integrable systems, the 
real parts of resonances, arranged in increasing order, 
form an essentially uncorrelated sequence. At fixed quan- 
tum number m, they form an ordered ladder, but the su- 
perposition of a large number of different, uncorrelated 
ladders results in a Poisson-like statistics. On the same 
grounds we assume uncorrelated resonance widths T, too. 
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The transmission amplitude at energy E from the 
lower m— channel to the upper I— channel is: 



S, n -,l+(E) 



k 2 (E,m)h 2 



a 2 e 2ik(E,m)L y 



(k(E, m)U + iaY x 



J ml i 



(4) 



where k(E, m) = \/2EhT 2 — 4ir 2 m 2 L x 2 . A computa- 
tion shows that the residue of (Q) at a resonance pole is 
~ — iT/2 at small T. Hence, the resonant part of the 
total transmission coefficient: 



T{E) = \Sn 



,771+ 



(E)\ 2 



(5) 



(the sum being over all open channels at energy E) has 
the form (|l|) apart from a slowly varying factor. The 
smooth part is oc E 5 ^ 2 /a 4 . Collecting various estimates 
we see that at h/(aL y ) < 1, h 2 / ' L 2 < FT < E < a 2 , 

E^ 2 L x h~ 1 a~ 2 3> 1 the assumptions of our general ar- 
gument are satisfied, hence the graph of T(E) vs E 
in (E - W < E < E + W) should be fractal with 
dimension / = 4/3, over scales intermediate between 

S mm ~ 1.8e 3/4 r^ o 4 and 6 max ~ r_E . 

A numerically computed graph of T(E) is shown in 
Fig.l (lower). The corresponding fractal analysis is 
shown in Fig. 2 and fully confirms the theory. In that 
case, the above estimates for 6 m i n , 5 max give 10 and 
10 3 ' 9 respectively. 



A generalization of the above model allows for the in- 
vestigation of parametric fluctuations. It is obtained by 
using boundary conditions -0(0, y) — e lc ^ip(L x , y). This is 
equivalent to a particle moving on a cylinder, with axis 
in the y— direction, enclosing a magnetic flux 4>. This 
problem is still completely integrable. Replacing m by 
m — 0/(2-7r) throughout the equations derived at 4> = 
yields the corresponding theory. At hxed energy E, the 
total transmission fluctuates as <f> is varied, in the manner 
illustrated in Fig.l (upper). The theory of such fluctua- 
tions is completely parallel to the one we have described 
for fluctuations vs energy at fixed <f> = 0. Resonances de- 
pend on 4>, z„ )m = z n ^ m (4>), and the complex values of 4> 
solving the equations z n , m (<^) = E define resonance poles 
in the complex 0-plane. Omitting computational details, 
the real parts of such poles are distributed in [0, 2ir] with 
a mean spacing ~ n hy/2/E/L x . The distribut ion 

1/3 

of their widths behaves like at small F^, with a 

mean ~ L 2 TEe4>/{8ir 2 h 2 ). Our general discussion is thus 
valid for parametric fluctuations, too. In fact numerical 
data shown in Fig. 3 demonstrate the predicted fractal 
dimension 4/3. 
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FIG. 1. Total transmission versus flux cj> at fixed energy 
E — 5 x 10 10 (upper plot), and versus energy E at fixed flux 
= (lower plot), for a = 3.5 x 10° and a — 10 4 respec- 
tively. In both cases % = 1, L x = 2, L v = 0.4. The energy 
range in the lower plot is (E , E + W) with E — 5 x 10 7 , 
IT = 2.5 x 10 5 . Units are arbitrary. 
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FIG. 2. Fractal analysis of the graph of T(E) vs E shown 
in the lower fig |lj The meaning of S,N(6) is explained in 
the text. The straight line corresponds to fractal dimension 
/ = 4/3. The inset shows the autocorrelation 1 — C(6)/C(0) 
vs S. The straight line has slope 4/3. 

Both for energy-dependent and for parametric fluctu- 
ations we have computed autocorrelations of the fluctua- 
tion graphs. For the case of fluctuations vs energy, such 
correlations are defined by 

C(6) = (T( fl \E)TW(E + 5)) E 

where T^ 1 ' is obtained from T by subtracting a smooth, 
slowly varying part, and the average is taken over the 
scanned interval of energies. Treating T^ l \E) as a st a _ 
tionary stochastic process, one easily finds that such cor- 
relations behave at small S like C(0) — const, x <5 27 , where 
7 is the scaling exponent of the rms increment of T(E) 
over intervals of length 6: (\T(E + S) - T{E)\ 2 ) E oc 6 2 ~< 
(note that subtracting the smooth part does not alter 
fractional scaling). Generally speaking, the rms incre- 
ment is a quite different quantity from the average excur- 
sion which enters the definition of the fractal dimension. 
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With the present strong statistical properties the two 
quantities scale in the same way at small S: 2j = 2 — a, 
so C(S) ~ C(0) — const, x 5 4 / 3 , as confirmed by numer- 
ical data. However fractional scaling of correlations is 
not in general a sufficient condition for fractality. For 
instance, correlations in Fig. 2 and in Fig. 3 exhibit a 4/3 
scaling down to small 5 scales below the fractal range. 
This is because they are still determined by the statis- 
tics of narrow, non-overlapped individual peaks, which 
do not produce fractality any more. On the other hand, 
on increasing S correlation functions depart from the pre- 
dicted fractional behaviour already at values well within 
the fractal range, because a larger statistics is needed for 
increments recorded over a finite 5— grid to sample the 
distribution of excursions over the same grid. 



mands careful analysis. 

(3) uncorrelated energy spectra also occur with fully 
chaotic systems in the regime of strong quantum localiza- 
tion. A prototype system in this class is the kicked rotor, 
the quasi-energy spectrum of which has a Poisson-likc 
statistics in case of strong localization |l2| . For this class 
of systems, the T (differential) distribution behaves like 
1/r down to very small scales |L^ , |l4| ]. This explains re- 
cent findings || of parametric fractal fluctuations with di- 
mension 3/2, detected in the survival probability at fixed 
time, on varying a magnetic flux <f>. 

Support from MURST Research Project "Chaos and 
localization in classical and quantum mechanics" is grate- 
fully acknowledged. 
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FIG. 3. Fractal analysis of the flux-dependent fluctuations 
shown in the upper fig. 1 . The straight line corresponds to 
/ = 4/3 . The inset shows the correlation scaling exponent 
2 7 = 4/3. 

The work reported in this Letter hinges on the fully 
general fact, that fractality of scattering fluctuations is, 
first and foremost, a matter of complex level statistics. 
As in the case of real level statistics (closed systems) , in- 
tegrable systems have special properties in this respect, 
and we have in fact demonstrated fractal reaction curves 
for a completely integrable process. A few remarks are 
in order about the relevance of this finding: 

(1) In completely integrable scattering processes there 
is no mixing of flux between different channels. Like 
pseudo-randomness of energy levels, fractality comes of 
superimposing non-fractal fluctuation patterns from dif- 
ferent, uncorrelated channels. Therefore the nature of the 
reaction curves depends on the relative weight assigned 
to different channels. 

(2) The present result may also be relevant to quasi- 
integrable systems, which possess large stable compo- 
nents in their phase space. Such components may con- 
tribute a significant set of resonances, produced by tun- 
neling through invariant manifolds, with the statistical 
properties considered in this Letter. Their interplay with 
critical structures at the border of the stable regions de- 



ll] Y.V.Fyodorov and H.J. Sommers, J.Math.Phys. 38 
(1997)' 1918. 

[2] U.Smilansky, in the Proceedings of the 1989 Les Houches 
Summer School on Chaos and Quantum Physics, 
M.Giannoni, A. Voros and J. Zinn- Justin eds., North Hol- 
land, Amsterdam 1991. 

[3] R.A.Jalabert, in the Proceedings of the International 
School of Physics "Enrico Fermi", course CXLIII, 
G.Casati, I.Guarneri and U.Smilansky eds,. IOS Press 
(2000), and references therein. 

[4] T.Geisel, A.Zacherls and G. Radons, Phys. Rev. Lett. 59 
(1987), 2503 

[5] R.Ketzmerick, Phys. Rev. B54 (1996) 10841; L.Hufnagel, 
R.Ketzmerick and M.Weiss, Europhys. Lett., 54 (2001) 
703. 

[6] G.Casati, I.Guarneri and G.Maspero, Phys. Rev. Lett. 84 
(2000) 63. 

[7] H.Hegger, BHuckestein, KHecker, M.Janssen, A. 
Freimuth, G.Reckziegel and R.Tuzinski, Phys. Rev. Lett. 
77 (1996) 3885; A.P.Micolich, R.P.Taylor, R.Newbury, 
J.P.Bird, R.Wirtz, C.P.Dettmann, Y.Aoyagi and 
T.Sugano, J. Phys: Condens. Matter 10 (1998) 1339; A.S. 
Sachrajda, R. Ketzmerick, C. Gould, Y. Feng, P.J. Kelly, 
A. Delage, and Z. Wasilewski, Phys.Rev.Lett, 80 (1998) 
1948. 

[8] E.Louis and J.A.Verges, Phys.Rev. B61(2000) 13014. 
[9] G.Benenti, G.Casati, I.Guarneri and M.Terraneo, Phys. 

Rev. Lett., 87 (2001) 014101 
[10] Transport fluctuations in integrable cavities have been 

studied in detail, but their fractality doesn't appear to 

have been investigated as yet. See, e.g., ref. ||. 
[11] M.V.Berry and M. Tabor, Proc. R. Soc. London, Ser. A, 

349, (1976) 101. 
[12] M.Feingold, S.Fishman, D.R.Grempel and RPrange, 

Phys. Rev. B 31 (1985) 6852. 
[13] G.Casati, G.Maspero and D.L.Shepelyansky, Phys. Rev. 

Lett. 82, (1999), 524. 
[14] M.Titov and Y.Fyodorov, Phys.Rev. B61 (2000) 2444; 

M.Terraneo and I.Guarneri, Eur.Phys.J. B18 (2000), 303. 



4 



